Effect of three-body loss on itinerant ferromagnetism in an atomic Fermi gas 
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A recent experiment has provided the first evidence for itinerant ferromagnetism in an ultracold 
atomic gas of fermions with repulsive interactions. However, the gas in this regime is also subject to 
significant three-body loss. We adopt an extended Hertz-Millis theory to account for the effect of loss 
on the transition and on the ferromagnetic state. We find that the losses damp quantum fluctuations 
and thereby significantly increase the critical interaction strength needed to induce ferromagnetism. 
This effect may resolve a discrepancy between the experiment and previous theoretical predictions of 
the critical interaction strength. We further illuminate the impact of loss by studying the collective 
spin excitations in the ferromagnet. Even in the fully polarized state, where loss is completely 
suppressed, spin waves acquire a decay rate proportional to the three-body loss coefficient. 
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The Stoncr transition from a paramagnetic metal to a 
ferromagnet is one of the earliest known and seemingly 
simple examples of a quantum phase transition. Yet re- 
cent theoretical work [l|, Q has revealed a great deal of 
complexity and suggested that quantum fluctuations play 
a vital role in determining the behavior near to the quan- 
tum critical point. Specifically fluctuations drive the fer- 
romagnetic transition first order at low temperature and 
may lead to the formation of novel phases [l| . Whether 
these effects can explain puzzling experimental observa- 
tions in materials such as ZrZn2 and Sr3pUi207 Q, or 
whether coupling to phonons or other auxiliary degrees 
of freedom is involved remains an open question. Ul- 
tracold Fermi gases tuned by a Fcshbach resonance now 
offer experimentalists unprecedented control over many- 
body phenomena and so provide a concrete platform from 
which to answer such questions and enhance our under- 
standing of quantum critical phenomena in itinerant fer- 
romagnets 0, H, 

However, the atomic systems are not free of complica- 
tions of their own. A major obstacle to the formation of 
a ferromagnetic phase in an ultracold atomic gas is the 
loss of atoms due to three-body interactions Q • A recent 
and seminal experiment Q overcame this hurdle with a 
non-adiabatic Fcshbach field quench to festinate into the 
ferromagnetic state, and has now provided the first firm 
evidence of the formation of a ferromagnetic state in an 
atomic gas 0. The experiment has suggested that the 
ferromagnetic transition takes place at kpa rs 2.2, which 
is in stark contrast to the Quantum Monte Carlo predic- 
tion of kpa w 0.85 [l[ or the analytical prediction with 
quantum fluctuations of kpa m 1.05 @. Although both 
the trap geometry and non-adiabatic conditions can ar- 
tificially raise the required interaction strength the 
source of the discrepancy in the interaction strength has 
yet to be adequately understood. It is vital to resolve this 
inconsistency to definitively prove that the ferromagnetic 
phase was formed rather than some alternative strongly 
correlated phase [l(| . 

In this letter we investigate how the three-body loss 



interaction impacts on the Stoner transition and the be- 
havior of the system in the ferromagnetic phase. The 
loss can be viewed as an additional interaction that may 
renormalize the value of the repulsive interaction, as well 
as introduce fluctuations of its own that give rise to new 
phenomena. For example, it was recently argued that 
strong two or three body loss in a Bose system can give 
rise to effective hard core interactions [HI, [12, EH , and 
in particular and lead to dynamic formation of a Tonks- 
Girardeau gas [ill ]. In the itinerant fermion system, wc 
discover that three-body losses inhibit quantum fluctua- 
tions out of the fully polarized state which, since quan- 
tum fluctuations promote ferromagnetism [l], [||, raises 
the interaction strength required to stabilize the ferro- 
magnetic state. This mechanism provides strong motiva- 
tion to study the collective modes of the fully polarized 
phase. We find that spin waves acquire a finite life-time 
in presence of loss, even in the fully polarized phase where 
loss is completely suppressed in the pristine system. 

The renormalization of interactions due to atom loss 
can be studied by adapting the formalism developed in 
Ref. [|| that demonstrated how quantum fluctuations 
drive the ferromagnetic phase transition first order. Since 
the prediction of the interaction strength kpa = 1.054 
for the onset of ferromagnetism is backed up by robust 
Quantum Monte Carlo calculations [H, this formalism 
provides a solid foundation for the present study. 

In order to study the effect of loss using a frame- 
work suited for equilibrium systems, we adopt a general- 
ized linear response strategy in which we seek the modi- 
fied quasi-equilibrium state established by the small loss 
term. To abrogate the three-body loss we insert an arti- 
ficial single-body atom source term. Then following the 
prescription laid out in Ref. @ we integrate out quantum 
fluctuations and find that the three-body loss renormal- 
izes the interaction strength, a consequence of atom loss 
damping the quantum fluctuations. Once the magnetiza- 
tion is established we focus on those regions of the phase 
diagram where atom loss and therefore the source term 
is zero and so faithful to the experimental setup. Having 
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laid out the strategy we now present the formal calcula- 
tion including atom loss. 

Our goal is to calculate the quantum partition func- 
tion with a fermionic coherent state path integral Z = 
fT>(tli,ip)e~ 3 ^'^. The starting point is the action 
S[)j\ i/'], which describes a two component Fermi gas with 
pscudospin a G {T, 1} that is represented by fermionic 
fields ip and ijj 



S=l drdr 



Vv (dr + ek-M) Vv+.gi/'T'01'0iV'T 



(1) 



where the atoms have the single particle dispersion 
6k = fc 2 /2 and feel a repulsive s-wave contact interac- 
tion gS 3 (r) that can be tuned via a Feshbach resonance, 
and f3 = 1/fcsT is the inverse temperature. We have 
also set h = m = 1. Before proceeding we consider 
how to incorporate atom loss into the model. If we mo- 
mentarily consider the non-interacting system, and in- 
sert a hypothetical imaginary term ixipcripa into the La- 
grangian, then the corresponding retarded Green's func- 
tion is Gj^ = (d T + £k — fi + i%) -1 . If we now perform 
the analytical continuation r i— > it into real time, and 
consider the Fourier transform to the temporal domain 
we recover G^(t) = exp[ii(ek — p) — xi]9(£). Therefore, 
within linear response theory, we can identify an imagi- 
nary term in the action with the effective lifetime x" 1 °f 
the state |k), and therefore a probe of atom loss. 

In the cold atom gas the dominant loss process is 
three-body interactions. Following the linear response 
template we incorporate three-body loss into the action 
through the term iApV'tV'j.V'j.V'T- Here we replaced a pair 
of dynamic fermion fields with the average density p, 
which amounts to a two-body mean-field approximation. 
The three-body loss term is now on equal footing to the 
two-body interaction, and so is perfectly poised to study 
the renormalization of the interaction strength. The co- 
efficient A is estimated, in the weakly interacting regime 
fcpa << 1, to be A = llle(fcFa) 6 0- Since we seek a 
quasi-cquilibrium state, we counterpoise the atom loss 
with an atom source term — vy ippipa-, where 7 will be 
determined later. Within this formalism, the action is 



S 



drdr 



^ Vv (d T + e k - p - 17) 



+ (.9 + iAp) V'tV'lV'lV'T 
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To proceed we calculate the free energy following the 
prescription laid out in Ref. @. First we introduce a 
Hubbard-Stratonovich transformation in both the den- 
sity channel p and the magnetization channel <f> to de- 
couple the quartic terms in the fermionic field. This 
leads us to identify the spectrum £k,o- = £k + 17k + (.<? + 
i\p)(p — <j(j)) — p. After integrating out the fermionic vari- 
ables we expand the fluctuations in the bosonic fields to 
quadratic order and also integrate them out. To remove 



the unphysical ultraviolet divergence of the contact inter- 
action we employ the standard regularization setting q 1— » 

~^r~v\~^ri z^k3,4^ k i.T+4k 2 ,i~4k3.T-4k 4 a) lm- 

The prime indicates that the summation is subject to the 
momentum conservation ki+k2 = k3+k4, the Fermi dis- 
tribution n is calculated in the presence of the now imag- 
inary chemical potential, and v is the density of states 
at the Fermi surface of an equivalent non- interacting gas. 
This regularization allows us to characterize the strength 
of the interaction through the dimensionless parameter 
k-pa, where fcp denotes the Fermi wave vector and a is the 
s-wave scattering length. The analysis yields a perturba- 
tion expansion in terms of the dimensionless interaction 
strength and the loss parameter 



S = X]( ek + i 7k)"(Ck, CT ) + 



2kpa 
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T. 



(3) 



where the quantum fluctuations are encoded in the term 
T= y^ XCkLtM^q) N£k 3 ,t) + w(£k 3 „j)] ,q 

t &i,T +£k 2 ,; - £k 3 ,t - £k4,J. 



A similar ex pres sion was derived in the homogeneous 
case in Ref. [15( using second order perturbation the- 
ory. Consistent with the initial identification of atom 
loss through linear response theory, and the perturba- 
tion expansion in the interaction strength, we now ex- 
pand the action out to quadratic order in 7 and A, 
and assume that the gas is at low temperature so that 
n(£ + A) « n(£) - A6(£) - A 2 <S'(£)/2. To establish the 
connection to the free energy we first demand that there 
is no imaginary component to the free energy which fixes 
the atom source term to 

Here v a denotes the density of states at the Fermi energy 
of species a. We note that the atom loss is zero when the 
system is fully polarized. The real part of the action 
corresponds to the free energy F = Fo + A(A) modified 
by atom loss. The standard theory of a lossless system 
is encoded in the term 



z? ic \ , 2kFa 1 ( 2k F a \ 
F Q = 2^ ekn(£k, CT ) + p^pi - 2 | | T 
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The renormalization of the interaction strength due to 
atom loss enters through 



A = 2A 2 p 2 (T 

+ ^P7 [P](y i - 2p, i v[) + pj_(i/ T - 2p ] v\)] 



v\ + v\, 



P\ v \ - H v i 



(7) 
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FIG. 1: (Color online) The required interaction strength to 
reach the shaded full polarization with changing loss rate A in 
(a) the mean-field case and (b) when fluctuation corrections 
are considered. The atom loss variation with weak interac- 
tions is shown by the red dotted line, and experimental [l6|] 
estimates of atom loss by the red points with error bars. In 
(a) the regions where collective modes are underdamped and 
overdamped are highlighted, in (b) the fully polarized bound- 
ary with 7 > (solid) and 7 = (dotted) is plotted. 



where v' a is the differential of the density of states at 
the Fermi surface of species a. In our further analysis of 
the free energy it will prove convenient to consider the 
dimcnsionless loss parameter A = -KvXp/2. 

We now use the formalism to determine the phase di- 
agram, focusing on the boundary of the fully polarized 
state. The mean-field case presented in Fig. [Ua) pro- 
vides a solid foundation on which to build the analysis. 
Quantum fluctuations destroy the ferromagnetic state, 
but atom loss damps fluctuations and stabilizes the fully 
polarized state, which consequentially can be seen at 
weaker interaction strengths. This behavior stems di- 
rectly from Eq. (JT]), which has a negative gradient with 
magnetization. However, Refs. [l|, 0] demonstrated that 
to higher order in the interaction strength quantum fluc- 
tuations in fact stabilize the ferromagnetic phase. There- 
fore, when higher order corrections are taken into account 
in Fig. [ljb), quantum fluctuations at A = reduce the 
critical interaction strength compared to the mean-field 
case. These quantum fluctuations also drive the tran- 
sition first order. With increasing three-body loss the 
quantum fluctuations are damped, which in turn raises 
the interaction strength required to sec fcrromagnctism. 



This feature is repeated in Eq. (JT]) which has a positive 
gradient with magnetization. Beyond a critical value of 
the loss coefficient A > 1.7 which is of the order of the in- 
teraction strength g = ?t/2 at the Stoner criterion for the 
second order transition, the fluctuations that drove the 
ferromagnetic transition first order arc sufficiently sup- 
pressed to recover that second order transition. At high 
interaction strengths the border adopts a linear behavior 
as the highest order terms in the free energy in inter- 
action strength and loss are both quadratic. However, 
in this regime it is likely that higher order corrections 
will play an important role. The curve was calculated 
with the atom source bath in place; we also verified the 
robustness of our result by making the alternative ap- 
proximation of explicitly setting 7 = and ignoring the 
atom loss. In this case the curve (also shown in Fig. [TJs) 
is found to be qualitatively the same, though with an 
enhanced critical interaction strength. 

We now turn to consider the experimental repercus- 
sions of the renormalization of the interaction strength. 
The theoretical prediction Q for the variation of the 
loss coefficient as the system is tuned closer to the Fes- 
hbach resonance is given by A = llle(fcFa) 6 . As seen 
in Fig. [ljb), this trajectory does not enter the fully po- 
larized regime, so according to this model the system 
remains steadfastly in the paramagnetic (or at most par- 
tially polarized) phase. However this prediction is only 
valid for kpa <C 1. Indeed the loss coefficient as esti- 
mated from the experiment (points with error bars) 
is much lower for given It-pa. With this estimate the sys- 
tem in fact passes into the fully polarized state. Though 
the error bars on the experimental results are large, the 
results indicate that the system becomes fully polarized 
at kpa ~ 1.9, which compares favorably with the ex- 
perimental observation that the atomic gas became fully 
polarized at kpa « 2.2 Q. The experimental observa- 
tion of the raised critical interaction strength is strongly 
suggestive of the important role that quantum fluctu- 
ations have to play. The sole focus on the fully po- 
larized gas will make useful predictions for experiment 
since the analysis on the homogeneous system showed 
that the system is partially polarized over the narrow 
range 1.05 < kpa < 1.11 [6j, compared with the broad 
ambit of interaction strengths that exist within the trap 
< kpa < 8 [3]. Therefore, the partially polarized regime 
not considered here fills only a thin shell within the trap. 

Having exposed the significant renormalization of in- 
teraction strength due to three-body interactions we now 
demonstrate that the dispersion of the collective modes 
is a cogent probe of the impact of loss driven damping 
on the quantum fluctuations. 

The dispersion of the collective modes can be studied 
by searching for the poles in the propagators. We search 
only for the collective modes in the fully ferromagnetic 
state where the longitudinal mode is gapped so we fo- 
cus on the transverse modes (0, q) that have the inverse 
propagator 1 + G T (w + fi, p + q)G l (w, p) fl. 

Equating this to zero yields collective mode frequencies 
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with both a real and an imaginary part. The real part, 
which gives the collective mode dispersion is 
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2 5 / 3 3 
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5fc F a 1 + A 2 /(A:Fa) s 



(8) 



Disregarding losses, the collective modes dispersion is 
identical to that found by Callaway 17|. The quadratic 
spin dispersion emerges as ferromagnetism breaks time- 
reversal symmetry, and is equal to that of a single mi- 
nority spin species particle propagating through a sea of 
majority spin particles. The dispersion rises with increas- 
ing interaction strength k-pa as the system becomes stiffcr 
against spin rotation. Atom loss introduces an additional 
energy penalty for fluctuations, consequentially the dis- 
persion also rises with the loss rate parameter A. To fully 
expose the influence that atom loss has over the disper- 
sion it is useful to focus on the instability to a partially 
polarized phase which develops at 



2 2 / 3 3 



2V318 
25 



A 2 



(9) 



Without three-body loss, the fully polarized phase be- 
comes unstable at kpa — 2 5 / 3 3/5 in accordance with the 
prediction of the mean-field Stoner model. At mean- 
field level quantum fluctuations destroy the ferromag- 
netic state, so Eq. ((9]) matches the boundary Fig. QJa) 
which demonstrates how increased loss reduces the re- 
quired interaction strength, whereas Fig.QTb) highlighted 
the opposite effect when quantum fluctuations are taken 
into account. Working at the mean-field level, there is a 
maximum loss rate, A = 2 2 / 3 3/5, beyond which the fully 
polarized state cannot be formed. 

In addition to renormalizing the dispersion, the pres- 
ence of a loss interaction also leads to the decay of the 
spin excitations. As the spin wave propagates the atom 
spins develop a component in the opposite spin direc- 
tion incurring atom loss, which in turn damps the spin 
wave. The characteristic inverse timcscale of damping, 
or width, of a transverse mode can be found from the 
imaginary component of its frequency, 
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q 2 2 5 / 3 3A 
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2 5(fc F a) 2 1 + A 2 /(A; 
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(10) 



We see that due to atom loss the spin- waves become res- 
onances that are characterized by a momentum indepen- 



dent quality factor 



^ ~ r 25/33 
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In Fig. HJa) we highlight the region Q < 1, where spin 
excitations completely lose their integrity. 

In experiment, these collective modes can be ex- 
cited and probed by spin-dependent Bragg spectroscopy. 
A variable wavelength optical lattice potential couples 
asymmetrically to the spin degrees of freedom, and 
thereby excites transverse magnetic fluctuations. The 
collective mode response could be studied through dy- 
namical fluctuations of the cloud spatial distribution as 
a function of wavelength, laser amplitude, and detuning. 

An experimental handle that could modify the atom 
loss rate would gift investigators the ability to fully ex- 
plore the consequences of atom loss. Though the atom 
loss rate cannot be reduced, it sets a base level that can 
then be artificially raised. An additional bosonic [HI or 
different fermion [19[ species could act as the third body 
in the atom loss process. The loss rate will be propor- 
tional to the density of this third species, which can be 
conveniently controlled. The latter case has already been 
investigated with the same lithium species employed in 
the pioneering ferromagnetism experiment [1, [ljj| . 

In this letter we have shown that three-body loss 
damps quantum fluctuations. This hinders the transition 
into the ferromagnetic state and is consistent with the 
experimental findings. Furthermore, we have highlighted 
signatures of this mechanism in the collective mode spec- 
trum. A novel phenomenology was developed to probe 
the consequences of atom loss, and its generality opens 
the possibility to explore transitions from a phase with 
atom loss to one without. Candidate systems include the 
BEC-BCS crossover, p-wave superfluids, and the boson 
atom-molecule superfluid transition [20] • This system 
has no base atom loss, but adding a third species will 
activate molecule loss from the atomic side of the transi- 
tion. Thus the loss rate can be tuned from zero upwards 
and the impact on the phase transition fully exposed. 
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